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We propose a new framework to represent the perturbative S-matrix which is well-defined for all 
quantum field theories of massless particles, constructed from tree-level amplitudes and integrable 
term-by-term. This representation is derived from the Feynman expansion through a series of partial 
fraction identities, discarding terms that vanish upon integration. Loop integrands are expressed 
in terms of “Q-cuts” that involve both off-shell and on-shell loop-momenta, defined with a precise 
contour prescription that can be evaluated by ordinary methods. This framework implies recent 
results found in the scattering equation formalism at one-loop, and it has a natural extension to all 
orders—even non-planar theories without well-defined forward limits or good ultraviolet behavior. 


INTRODUCTION 


Since the revolutionary developments in quantum field 
theory around the middle of the last century, Feynman 
diagrams have been an essential tool for computing scat¬ 
tering amplitudes. As represented by the Feynman ex¬ 
pansion, amplitudes can be determined perturbatively by 
summing over all graphs with a fixed number of loops 
that connect the external states by appropriate vertices 
and propagators and integrating over all the internal 
loop-momenta. While intuitive and ultimately correct, 
the Feynman expansion rapidly becomes intractable both 
because the number of diagrams grows rapidly with the 
number of external legs and loops, but also because indi¬ 
vidual diagrams in many theories depend on a great deal 
of purely theoretical data (such as gauge redundancies) 
that have no consequences for physical predictions. 

Because of this, there has long been interest in finding 
formulations of perturbation theory without any explicit 
reference to unphysical data. This was the principal mo¬ 
tivation behind the Feynman tree theorem [T] , the essence 
of which can be be understood (in a rather novel way) 
from the well-known partial fraction identity: 
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Di - ■ ■ Dm 



( 1 ) 


When the factors I/Di are Feynman propagators for an 
off-shell loop momentum £, the terms (Dj-Di) become 
linear in £ and can therefore be interpreted as propaga¬ 
tors involving an on-shell momentum £. Expanding every 
term in the Feynman expansion in this way, it is tempt¬ 
ing to view the coefficient of each 1/Di{£) as a lower- 
loop amplitude, evaluated in the forward-limit. This is 
precisely correct for theories (such as those with super- 
symmetry) whose amplitudes are finite in the forward- 
limit; and this leads, for example, to the representation 
found in ref. [2] for one-loop amplitudes in supergrav¬ 
ity, derived within the scattering equation formalism [3]. 


But for most theories, however, scattering amplitudes di¬ 
verge in the forward-limit, preventing any (unqualified) 
interpretation of the terms in this expansion as limits of 
lower-loop amplitudes. 

In this Letter, we present a new representation of am¬ 
plitudes which avoids this obstruction, allowing one to 
write the coefficient of each off-shell propagator 1/Di{£) 
in the expansion above in terms of “Q-cuts”, which are 
well-defined in any theory in terms of gauge-invariant 
tree-amplitudes alone. In the following section, we de¬ 
rive this representation in detail for one-loop amplitudes 
and show how each term in this expansion can be inte¬ 
grated by ordinary means. We then describe how this 
approach generalizes to higher loops, illustrating the rich 
structure that emerges beyond two-loops. 

THE Q-CUT REPRESENTATION AT ONE-LOOP 

For the sake of concreteness and clarity, let us restrict 
our attention to one-loop amplitudes of theories with 
massless particles, so that Feynman propagators for the 
loop momenta are of the form N{£)/{£ +P)'^. Both in 
order to generalize the partial fraction decomposition to 
cases with loop-dependent numerators, and to introduce 
ideas that will prove useful later on, let us derive equation 
([^ as an instance of Cauchy’s residue theorem. Consider 
transforming £r^£ + rj where rj is orthogonal to the exter¬ 
nal momentcQ satisfying rf = z so that + z, and 

Di^^ Di + z for propagators involving £. Transforming a 
Feynman diagram this way and dividing by z results in 
N{£) 1 N{£,z) 

D^{£y ■ ■ Dm{£) ^ ~z{Di{£)+z)---{Dm{£)+zy 
After this deformation, the statement that the sum of all 
the residues in z vanishes trivially reduces to equation 
0 as a special case. 


^ This is always possible in the framework of dimensional regular¬ 
ization by taking rj to lie along the extra dimensions. 
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Using this generalization of the partial fraction expan¬ 
sion to decompose every integral of the Feynman expan¬ 
sion, it is easy to see that every term has the form (up 
to a shift in i by external momenta), 
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(2UPi + Qi)---(2UP„-t-g„) 


(3) 


where iV(£) accounts for both the numerators of the di¬ 
agrams and any loop-independent propagators. Let I(^) 
denote the factors in the square brackets above. 

Although forward-limit divergences prevent us from in¬ 
terpreting I(^) in terms of an entire tree-amplitude in 
general, it turns out that we can construct T-it) in terms 
of tree-level objects up to terms that vanish upon inte¬ 
gration. This becomes possible after one further partial- 
fraction-like expansion—this time, in the scale of i. Con¬ 
cretely, consider the residue theorem resulting from: 


represent graphically as follows: 


(7) 


Notice that the shifted propagator, corresponding to the 
factor 1 /{{i + Pl)'^ -i'^) in equation (|^, is indicated by a 
dashed line in the figure above to distinguish it from the 
unshifted, off-shell propagator l/£^. 

We claim that the sum over all Q-cuts (with 7 ^ 0) 
reproduces any one-loop amplitude. Notice that the inte¬ 
grand of a Q-cut is similar to a Cutkosky unitarity cut [1] . 
The principal novelty involved in the Q-cut is that the 
amplitudes involved are evaluated with shifted (on-shell) 
values of multiplied by unusual propagators. 



I(£)^I(a,£)^^^. (4) 

Clearly, I{i) is recovered as the residue of J(q:, £) at a = 1; 
and by Cauchy’s theorem, this is equal to (minus) the 
sum of all other residues. These residues are associated 
with three types of poles: at zero, at infinity, and at finite 
locations (a^f 1). By inspection of equation ([^, residues 
at q: = 0 correspond to integrals of the form, 

k 

where the product runs over only those factors for which 
Qk = 0 , and N{£) denotes all other factors at a = 0 . 
Integrals of this form must vanish upon integration in 
any number of dimensions d because the denominator is 
homogeneous in i and hence scale-free. Similarly, the 
Laurent expansion of I{a,i) at a —>00 can involve only 
terms homogeneous in i, which hence vanish upon inte¬ 
gration. Notice that these residues precisely correspond 
to the terms poorly defined in the forward-limit. 

Therefore, we can replace X{£) by the sum of residues 
ofX(Q;,^) at a^{0, 1 , 00 }. Importantly, all such residues 
can be interpreted as involving two additional on-shell 
particles, with specific momenta determined by the suc¬ 
cessive residues. It is not hard to see that expanding 
every term in the Feynman expansion in this way, the 
coefficient of each pair of propagators becomes a prod¬ 
uct of complete tree-amplitudes evaluated for particular 
on-shell, internal momenta (and summing over states): 

Jl, -£r)^ (2£-Pl-\- Pi) ( 6 ) 

with £l = ol{£ rj) and £r = £l+Pl, with rf = —£'^ 
and a= - Pi/{2£-Pif)^Q, and where Pl denotes the sum 
of momenta over a partition of external legs. We refer 
to functions of the form of as Q-cuts, which we can 


CONTOURS OF INTEGRATION 

At a fundamental level, the causal structure of 
scattering amplitudes is encoded in the Feynman ie- 
prescription, critical to the precise definition of the loop 
integration contour. It will be useful here to observe 
that every Feynman propagator can be assigned its own 
e, transforming Dj 1 —^ Dj +iej] so long as each Cj is real 
and positive, the physical contour will be unambiguous 
(and independent of the Cj’s). 

Since Q-cuts do not involve products of Feynman prop¬ 
agators, it is not immediately clear how to assign fe’s 
to the linear poles appearing in their definition. How¬ 
ever, if we had started with a single Feynman integral 
with specific e’s for each propagator, then the partial 
fraction expansion would result in terms of the form 
Dij = (Di-Dj), with contours prescribed by shifts in¬ 
volving icij = i{ei-ej), the signs of which will be fully 
determined by the (arbitrary) ordering of the original 
e’s. This always provides a precise contour of integration 
for the resulting expressions that is guaranteed to match 
the original expression. The problem is in going in the 
other direction: to assign an unambiguous prescription 
for the e’s associated with the linear-factors in £ of each 
Q-cut integral}^ 

(We should mention that the contour we describe here 
requires that on-shell tree amplitudes are represented in 
a way that involves only local poles. Representations of 
trees generated by the BCFW recursion relations [S], for 
example, involve spurious, complex poles in individual 
terms. Finding a contour prescription for such terms is 
an important and interesting open problem.) 


^ In planar theories, where all the e’s of all propagators can be 
identified and conventionally ordered according to the ordering 
of the external legs, an unambiguous convention for the signs of 
the e’s for any linear factor can be easily assigned. 
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It turns out that a maximally democratic contour pre¬ 
scription will always work—meaning, that it is guaran¬ 
teed to match the Feynman contour. Specifically, we may 
simply average over the possible relative orderings of the 
Cj’s for the initial propagators. 

It will help to clarify this contour prescription with a 
few examples. Consider a single propagator linear in €, 
denoted Xjx. Its e can be arbitrary, so we are instructed 
average over the 2! possible choices for the sign of its le- 
prescription; thus, the contour prescription becomes the 
principle value via the replacement: 


1 1 

- 1 -^ - 
X 2 
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x + if. x-ie 


= V 


LxJ 


( 8 ) 


(Recall the possible contour prescriptions—defined by. 
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a; ± je 


= V 
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=F iirS^x), 


(9) 


where S(x) is the (non-holomorphic) Dirac d-function.) 

For three propagators (the case involving two linear 
propagators after the partial fraction decomposition), the 
result is rather less trivial: the 3! relative orderings for 
the e’s of the three propagators result in a more inter¬ 
esting distribution of differences. It is a simple exercise 
to see that the resulting contour prescription for the any 
pair of differences should be: 
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xy 

.X. 


-y- 


- ^7r2(5(x)(5(2/). 


( 10 ) 


For four propagators (three linear ones), the prescription 
is quite similar: l/(xyz) should be replaced by. 
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1- 

V 
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V 
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( 11 ) 
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(S(x)6(y)V - +6(y)6(z)'P — +S(z)6(x)V 
V Izi La:J 


1 

LyJ 


And for five or more, there will be terms involving four 
^-functions for each quadruplet (times -htt'^/S), etc. 


ONE-LOOP EXAMPLES 


To exemplify the contour prescription, consider a sim¬ 
ple (massive) bubble integral with two propagators. Af¬ 
ter partial fractioning, and using the contour prescription 
described above, it becomes: 


( 12 ) 


To simplify the ^integration, it is possible to combine the 
denominators via Schwinger parameters, as usual. The 
following identities, which generalize the usual distribu¬ 
tion formula, i/{x + ie)^ are useful: 

r2?i r°° 

da sign(a) 2ttS{x) da (13) 


rd‘^e 

1 ^ 

2 

j T^d/2 

£2 -b it 

.2£-p -bp2. 


V 




Integrating over £ and an overall Schwinger parameter, 


the integral (121 readily becomes: 

sign(a) 


F 2- 


da 


(—a(l — a)p'^ — ze)2“‘^/2 


(14) 


One can see that the regions a > 1 and a < 0 exactly cancel 
each other, leaving the region 0 < o < 1, which precisely 
reproduces the usual Feynman-parameter expression for 
this integral. We could similarly prove, by exploiting can¬ 
celations in Schwinger parameter space, that any Feyn¬ 
man integral is correctly reproduced by the sum of its 
Q-cuts, integrated using the contour prescribed above. 

As a relatively simple but illustrative application, con¬ 
sider the amplitude for incoming gluons with the same 
(+)-helicity in planar Yang-Mills theory. Considering for 
simplicity the contribution from a complex scalar loop, 
and substituting the product of trees using the notation 
of ref. [6], the Q-cut defined in equation ® becomes: 


[ 12 ] p^-£^\ 1 [pI,/{2£-p,,)Y 

(12) (2Api)J £2 (2Api2+p2^) 


[34] 


,(34) {-2£-pi)J ■ 

(The full amplitude will be represented by the sum of this 
Q-cut and its three cyclic rotations.) For this helicity 
amplitude, an identical result would be obtained for a 
gluon in the loopj^ Up to integrals that vanish upon 
integration, the preceding is equivalent to. 


[12] [34] 




(12)(34) £2(2£.pi)(2fpi2+p22)(-2£.p4) ' 

In the conventional unitarity method [3 i, the two- 
particle cut of this is recognized as that of a box integral 
with numerator, which is then integrated. Similarly, 
its single-cut matches that appearing in the Feynman 
tree theorem (see e.g. iiiin])- The distinctive feature of 
the present approach, however, is that each term can be 
integrated directly, bypassing the reconstruction of the 
off-shell integrand, and without requiring a forward-limit 
interpretation of the single-cut. 

Specifically, paying attention to the contour in equa¬ 
tion 0 and applying Schwinger parameters ( [T^ to the 
linear propagators, the integral over the Q-cut gives: 

-i [12] [34] |■°° dadbdc[l + Oid-4)] 

1^(12) (34) 

1 


{b{l-a-b-c)s + act — ie)2-d/2 


g l^sign(a)sign(6)sign(c) -b ^ (sign(a) -r sign(6) + sign(c)) 

where s = Pi2 and t = pf3 are the usual Mandelstam in¬ 
variants. After rescaling b by jcj, the a and c integrations 
both become elementary, and the b integration can be 
done in the limit of d—>-4. With this, we find the all-plus 
helicity amplitude to be given by: 

-i [12][34] / s . St 


12-167r2 (12) (34) \s + t ^ {s + \ t) ) 

Notice that logarithm nicely cancels in the sum, leaving 
the correct answer (see e.g. [6]), 
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Ai-i°°P( + ,+,+,+)=-- 


[12] [34] 


6 167r2 (12) (34) ’ 


(17) 


^ So long as only physical (d-dimensional) gluons are included in 
the polarization sums, then no ghosts are needed. 

















































4 


Regarding extensions to higher multiplicity and higher 
loops, it is worth mentioning that the Q-cut representa¬ 
tion can be combined with other modern techniques—for 
example, integral reduction and the use of integration-by¬ 
parts identities. 


EXTENSIONS TO HIGHER LOOPS 

To generalize the construction to two loops, we be¬ 
gin by writing each Feynman diagram such that only 
loop momenta ^1,^2 or (£1+^2) enter propagators. We 
then separately partial-fraction-out the propagators of 
each three type. More precisely, we introduce a three- 
parameter deformation rji where r]l = zi, 77! = Z2 

and (771 +772)^ =-Z3- Partial-fractioning then expresses the 
amplitude in terms of its residues in zi,Z2, Z3. 

In each variable Zi there are residues at both finite and 
infinite locations. The residues with all three Zi finite are 
immediately given by on-shell three-particle cuts, given 
by the following Q-cut: 



L 


R )<^ 




( 18 ) 


where £3 = ^1-1-^2 + ^*7: and the extra-dimensional com¬ 
ponents ti = di + rji are such that = 0 for i = 1 , 2 , 3 . 
The residues at infinity require more work. It is easy to 
see that residues with two or three Zi at infinity yield 
vanishing integrals and so can be discarded. Residues 
with one at infinity, say 23, represent degenerate topolo¬ 
gies with disconnected loops. For these terms we repeat 
the procedure introduced already at one-loop and rescale 
ii —>■ aiii and expand by partial fractions separately in 
ai,a2- This yields a second type of Q-cut: 



These tree amplitudes are evaluated using £i = ai{ii + rii), 
with rjf = and ai_2= - PI^h /similar to 
the one-loop case. The amplitudes are projected onto the 
Oth-order term in the Laurent expansion for large 771-772. 

Intuitively the two Q-cuts above account for graphs 
with connected and disconnected loops, with the extra- 
dimensional deformation and 771 • 772 —> 00 projection ac¬ 
complishing a gauge-invariant separation between them. 
According to our derivation, adding all Q-cuts (for all 
possible external leg insertions) will reproduce the correct 
integrated amplitude. The contour for each Q-cut is de¬ 
termined by equations ([^-([Tl]), which should be applied 
separately for each group of linear propagators—that is, 
those involving (A-ii+B), {A-i2 + B) or {A-{£i+£2) + B). 

As a natural variation of the same technique, in a pla¬ 
nar theory one could restrict to a two-parameter defor¬ 
mation with 23 = 0, since no graph has more than one 


mixed propagator l/{£i+£2 + P)^. The residues would 
then be related to the double-forward limits of trees, in 
theories where this limit can be defined m- 

We believe that the Q-cut-construction generalizes 
straightforwardly to any loop order. For example, using 
a 6 parameter deformation with all rji-rjj taken indepen¬ 
dent at three loops, we obtain the six Q-cut diagrams: 



SUMMARY AND CONCLUSION 

The desire to represent loop amplitudes directly in 
terms of lower-loop amplitudes is motivated by the prac¬ 
tical and conceptual advantages of eliminating explicit 
reference to the redundancies required by the Feynman- 
diagrammatic expansion. Recently, several such repre¬ 
sentations of loop amplitudes in terms of trees have ap¬ 
peared in the context of the scattering equation formal¬ 
ism. This is in part because this formalism makes it pos¬ 
sible (at least for certain theories) to systematically reg¬ 
ulate the divergences of tree amplitudes in the forward- 
limit [niisig But it remains an important, open prob¬ 
lem to systematically regulate the forward-limit diver¬ 
gences of amplitudes in general theories. 

In this Letter we have described a new, “Q-cut” repre¬ 
sentation of loop amplitudes, derived from general field 
theory arguments and without any reference to forward 
limits. And this representation naturally extends to all 
orders of perturbation theory opening new possibilities 
for computation. 

We thank Zvi Bern, Lance Dixon, Peter Goddard and 
Lionel Mason for helpful discussions. This work has been 
supported in part by a MOBILEX grant from the Dan¬ 
ish Council for Independent Research (JLB) and by Qiu- 
Shi funding and Chinese NSF funding under contracts 
No. 11031005 , No .11135006 and No.lI 125523 (BF). BF is 
grateful for the hospitality of the NBIA. 


[1] R. P. Feynman, Acta Phys. Polon. 24, 697 (1963). 

[2] Y. Geyer, L. Mason, R. Monteiro, and P. Tourkine, 
(2015), arXiv:1507.00321 


^ And in the case of planar i/j^-theory, such a prescription exists to 
all orders B], making it possible to write all L-loop amplitudes 
in terms of L-fold nested forward-limits of trees. 













5 


[3] F. Cachazo, S. He, and E. Y. Yuan, Phys. Rev. Lett. 
113, 171601 (2014), arXiv: 1307.2199 

[4] R. E. Cutkosky, J. Math. Phys. 1, 429 (1960). 

[5] R. Britto, F. Cachazo, B. Feng, and E. Witten, Phys. 
Rev. Lett. 94, 181602 (2005), arXiv:hep-th/0501052 

[6] Z. Bern and A. G. Morgan, Nucl. Phys. B467, 479 
(1996), arXiv:hep-ph/9511336 

[7] Z. Bern, L. J. Dixon, D. C. Dunbar, and D. A. Kosower, 
Nucl. Phys. B425, 217 (1994), arXiv:hep-ph/9403226 

[8] Z. Bern, L. J. Dixon, D. C. Dunbar, and D. A. Kosower, 


Nucl. Phys. B435, 59 (1995), arXiv:hep-ph/9409265 

[9] S. Catani, T. Gleisberg, F. Krauss, G. Rodrigo, and J.-C. 
Winter, JHEP 09, 065 (2008), arXiv:0804.3170 

[10] E. W. Nigel Glover and C. Williams, JHEP 12, 067 
(2008), arXiv:0810.2964, 

[11] S. Garon-Huot, JHEP 1105, 080 (2011). 

[12] G. Baadsgaard, N. E. J. Bjerrum-Bohr, J. L. Bourjaily, 
P. H. Damgaard, and B. Feng, (2015), arXiv:1508.03627 

[13] S. He and E. Y. Yuan, (2015), arXiv: 1508.06027 

[14] G. Baadsgaard et al, unpublished. 


